Anomalous Kondo Spin Splitting in Quantum Dots 
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The Zeeman splitting of localized electrons in a quantum dot in the Kondo regime is studied 
using a new slave-boson formulation. Our results show that the Kondo peak splitting depends on 
the gate potential applied to the quantum dot and on the topology of the system. A common fact 
of any geometry is that the differential susceptibility shows a strong non linear behavior. It was 
shown that there exist a critical field above which the Kondo resonance is splitted out. This critical 
field rapidly diminishes when the gate potential is lowered, as a consequence of the reduction of the 
Kondo temperature and a subsequent strong enhancement in the differential susceptibility occurring 
at low fields. The critical field is also strong depedent on the topology of the circuit. Above this 
critical field the Zeeman splitting depends linearly upon the magnetic field and does not extrapolate 
to zero at zero field. The magnitude of the Y-intercept coordinate depends on the gate potential but 
the slope of this function is not renormalized being independent of the value of the gate potential. 
Our results are in agreement with very recent experiments. 
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The observation of Kondo effect in quantum dots has 
already a large history 0, • The effect appears when a 
net spin becomes localized at a quantum dot (QD) cou- 
pled with metallic leads 

ill 

Local spin fluctuations 
(LSF) in the vicinity of the QD allows the spin of the 
electrons in the leads to screen the localized spin, as a 
magnetic impurity is screened in a host metal. A sharp 
resonance in the local density of states (LDOS) is devel- 
oped at the Fermi level of the system. When the QD 
is connected as a bridge between two metallic leads, the 
embedded geometry (EG) , the resonance provides a chan- 
nel of conduction and the conductance reaches the value 

2 

G = 2Go, being Go = t- the quantum of conductance. 
When the QD is side connected (SCG), two channels of 
conduction interfere destructively (one direct and other 
across the QD, from lead to lead) and the conductance is 
zero for the value of the gate potential, for which the in- 
terference is complete. The existence and characteristics 
of the Kondo resonance can be studied through conduc- 
tance G of the system. 

To better understand the Kondo effect it is important 
to study its behavior when an external magnetic field is 
applied, ft is expected the LSF to be quenched and the 
Kondo regime to disappear for large enough fields. For 
intermediate values, the Kondo resonance is splitted and 
its effect can be observed in the differential conductance 
of the system Q ■ In recent experiments 0, 0| it was pos- 
sible to obtain very precise value of A, the Zeeman split- 
ting, (ZS) for localized electrons in a QD in the Kondo 
regime (KR). However, these experiments are controver- 
sial. In one of them, 0, above a critical value of the mag- 
netic field, it was found that A ~ Ab = 2g^sB, where 



g, fiB and B is the gyromagnetic factor, the Bohr magne- 
ton and the magnetic field respectively. In others, U|(J, 
bellow a critical value, the effect of the magnetic field was 
not observed and above it the splitting was larger than 
Ab and does not extrapolate to zero at zero field. 

The splitting was initially theoretically predicted to 
be ±g/Af)B, in line with 2]. However, more recent calcu- 
lations .8] shows that the splitting would be observable 
only above B c , when the magnetic Zeeman energy As 
becomes competitive with the Kondo temperature Tk- 
Moore and Wen |i| established the relationship between 
the field induced spin splitting in the spectral function 
of the system in equilibrium and the splitting appear- 
ing in the differential conductance when the system is 
slightly out of equilibrium. They predict also that the 
spin screening characteristic of the KR should reduce 
the magnetic field splitting, i.e A < Ab- Unfortunately 
these theories were not capable of explaining the phe- 
nomenology seen in the latest experiments 0,u|. 

In this letter we developed a theoretical study that cor- 
rectly explains these last experimental results. We use 
different approaches to study the problem. We propose a 
new slave-boson formalism(SB) that considers two inde- 
pendent boson fields, each one associated to up and down 
spin. We solve the problem using as well the embedded 
cluster method (ECM), which is a powerfull tool to study 
strong correlated systems and we compare the re- 
sults of these two approaches. We study the system for 
the two geometries mentioned above (side connected and 
embedded) . The value of the splitting results to be de- 
pendent on the geometry of the circuit and has a strong 
dependence on the gate potential applied to the QD and 
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FIG. 1: Sketch of the geometry of the two topologies consider 
in this work. Upper panel, SCG. Lower panel, EG. 



the coupling between the leads and the QD. We also ob- 
tain the solution in the Hubbard I (HI) aproximation [lOj] 
that adequately describes the situation for T > T}.. 

The system is described by the Anderson impurity 
Hamiltonian including the effect of an external magnetic 
field. We identify three contribution to the total Hamil- 
tonian: the leads (Hl), the QD term (He) and the con- 
nections among the leads and the dot(-ffr)- Hl is the 
Hamiltonian of two semi-infinite chains for the EG con- 
figuration or of a one dimensional infinite chain for the 
SCG case. The Hamiltonian He is given by, 



H c = 



T n 0;CT + Un 0l1 n 0ll 



(1) 



where the subindex denotes the QD site, e±i = V g ± B 
is the spin dependent QD energy level, V g is the gate 
potential at the QD site and U the Coulomb repulsion 
that, for simplicity, is supposed to be infinite. We have 
adopted \gfis\ = 1- Finally, 



H T 



\Va\. a a 0o 



+ c.c. 



(2) 



where i is —1 and 1 for EG or i = 1 for the SCG. In the 
first case — 1 and 1 denotes the first site of the left and 
right semi-infinite one dimensional leads. In the second 
case 1 is the site of an infinite one dimensional lead to 
which the QD is connected and V is the hopping proba- 
bility(see Fig. 1). All the energies are expressed in units 
of the hopping among the leads sites t, D denotes the 
bandwidth of the leads and the Fermi energy is Ef = 0. 

In the new slave-boson formulation we define a boson 
field associate to each electron spin. The fermion opera- 
tors a\. a is written as the product, al. a = cJ. CT &t, where 
c and b denote the quasi-fcrmion and boson operators. 
Following a similar formulation to the standard one |l2| , 
Ht is transformed into 
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(3) 



In order to avoid double occupancy at the QD, constraint 
conditions have to be imposed. They are written, 
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(4) 



When both spins are equivalent, the above equations 
merge into the constraint of the standard slave-boson 
formalism. The Hamiltonian is treated in the mean-field 
approximation for the boson fields (MFSB), substituting 
in equations (3) and (4) the operator b a by its mean- 
value z a =< b a >. The mean- field versions of equations 
(4) are incorporated into the Hamiltonian using the La- 
grange multipliers "f a . Minimizing the total energy with 
respect to these multipliers we obtain two new equations, 



KV < a|. CT c 0CT > +zl(j a + 7- 



') +l<y -7- 



(5) 



where K = 1 for the SCG and K = 2 for the EG. We 
solve equations (5) self-consistently and obtain the Green 
function at the QD and the transmission T(E) across it 
(G = T(E f )). 

We have chosen the parameters to obtain a similar 
Kondo resonance for the SCG and the EG configura- 
tions. They are V = 0.35 for the EG and V = 0.65 
for the SCG. Increasing V g the system goes from Kondo 
to mixed valence regime. 

Fig. 2 (upper panels) shows the transmission at zero 
field obtained using the approximations MFSB and 
ECM. The effect of the Kondo resonance is clear: G = 2 
or G = for the EG or the SCG configurations respec- 
tively, and a resonance or anti-resonance width roughly 
similar for both geometries and both methods of calcu- 
lations. When the magnetic field is increased the Kondo 
peak is splitted and eventually the system goes out of 
resonance. In the medium panels of Fig. 2 we show the 
transmission for B — 0.18. The value of the magnetic 
field splitting obtained by the two methods are very sim- 
ilar, and for the ECM the results are almost independent 
of the cluster size exactly diagonalized. 

To study how this splitting develops we write down 
the expression for the Green function at the QD. For the 
MFSB solution we have, 



E 



2 la + K{Vz a fg L {E) 



(6) 



Here gh{E) is the undressed Green function of the lead. 
From equation JHJ we can obtain E a , the values of the 
renormalized energies of the QD level. They are, 



E„ 



1 - {{K - \)VzJtf 



(7) 



FIG. 2: Left column EG and right column SCG. Upper panels: 
transmission as a function of energy for B = 0. Continuous 
line SB, broken line ECM(11 sites). Medium panels: the same 
as above but for B = 0.18 and two different cluster sizes. 
Lower panels: the same as above for three different magnetic 
fields Bx < B 2 < B 3 obtained with MFSB. For more details 
see the text. 
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can be interpreted as the 



i-{K-l){Vz„f 

spin dependent renormalization of the resonance level by 
the interaction. Depending on the geometry (K = 1 or 
2) there is also a multiplicative correction in the term e CT . 



The ZS is given by the difference A = — E 



A = 



V g + B + 2 7i 



V g - B + 2 7T 



i-(K-i)(v Zl y i-(K-i)(vz r y 



(8) 



In Fig. 3 and 4 (lower panel) we show the parameters 
of the problem as a function of the magnetic field. For 
both configurations z\ — > 1 ,Z\ — > when B — > oo. A 
and Ar = Tj — Tj show a rapid increase at low fields, 
more evident for the SCG case. Ar results to be al- 
most constant for a large range of values of the magnetic 
field. As a consequence, A(B) tends asymptotically to a 
linear function. The slope of this linear function is qual- 
itatively the same as for the non interacting regime (for 
the EG configuration there is a small correction propor- 
tional to V 2 with respect to the non- interacting case). 
For low fields, A rises due to the Zeeman effect but also 
because the local level at the QD suffers a spin depen- 
dent renormalization as the magnetic field progressively 
quenches the LSF, giving an additional contribution to 
the splitting. The differential magnetic susceptibility 
(DMS) is strongly enhanced for low magnetic fields, (see 
the curve no,f — no,x = A n in figures 3,4). For larger 
fields, when the spin fluctuations are quasi totally sup- 
pressed, A recovers its linear function appearance, but 
shifted up with respect to the origin, due to the Coulomb 
interaction acting on each spin level. Lowering V g , the Y- 
intercept coordinate of this linear function is increased, 
although the slope does not change. For a fixed value 
of the field B, A results to depend linearly on V g with 
a slope a almost independent of B (see the inset of fig- 
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FIG. 3: Results for EG. Upper panel: The function A(B) for 
three different values of V g . The semi-sum of the resonance 
widths W, for Vg = —0.5, is shown with thick point line. The 
equation W(B) = A(B) defines B c . The point line shows A 
using the Hubbard I approximation. Lower panel: evolution 
of the parameters with B, for V g = —0.8. The inset shows 
A(V g ) for two magnetic fields. 



ures 3 and 4). This implies that, as known A in- 
creases logaritmically by lowering the Kondo tempera- 
ture since V g ~ {V 2 / D)log{D/T k ). As A > A B , the 
slope a is larger than the value predicted in and in 
agreement with recent experiments |3 16| . Increasing V g , 
we enter into the mixed valence regime and A recovers 
a nearly linear evolution that almost extrapolates to the 
origin. In this regime Ar is not constant anymore but 
increases monotonically with the field. As a consequence, 
the slope of the linear relationship between A and B is 
renormalized for both geometries with respect to the non- 
interacting case. In spite of this, it is important to realize 
that A is always lower than the value it acquires in the 
Kondo regime. 

To determine when the ZS can be experimentally de- 
tected the relevant quantity to be study is W, the semi- 
sum of the resonance width of each Kondo peak. It can 
be expressed as, 



W = K 



{Vzif 



\-{K-l){V Zi ) 2 l-{K-l){Vztf 



(9) 



The ZS is not detectable if A < W. The value of the 
critical magnetic field B c satisfies the equation A(B C ) = 
W{B C ). In figure 3 and 4 we plot W{B) for V g = -0.5 
and V g = —0.4, respectively. The values of W obtained 
from (9) at B = are the same to those obtained with 
the usual formula T k ^Dexp{-xV g /2KV 2 )^ where x is 
7T for EG or 27r for SGC. The B c obtained for each config- 
uration corresponds to a A(B C ) << T k in line with the 
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FIG. 4: Results for the SCG. Upper panel: The function 
A(B) for three different values of V a . The semi-sum of the 
resonance widths W , for Vg = —0.5, is shown with thick trace 
line. The equation W(B) — A(B) defines B c . The point line 
shows A using the Hubbard I approximation. Lower panel: 
evolution of the parameters with B, for V g = —0.5. The inset 
shows A(V g ) for two magnetic fields. 



experimental observation. Note that B c is appreciably 
smaller for the SCG than for the EG topology because 
in this case A has a negative contribution comming from 
V g (8), which is absent for the SCG, and also Ar has a 
more rapid increasing with B for the SCG. In figure 2 
(lower pannels) it is shown how the splitting develops for 
each geometry. 

We have also solved the problem in the Hubbard I 
approximation [ljj. It describes the situation for T > 
Tk, i.e. when the system is outside the Kondo regime. 
The local Green function results to be: 



1 



E-e a + KV 2 {l-n- a )g L (E) 



(10) 



We can determine the position of the resonances as a 
function of the external field llfj. A results to be: 



A 



V +B 



V a — B 



1 - (K - 1)V 2 (1 -n T ) 1 - (K - 1)V 2 (1 -ni) 



(11) 



The function A(B) is plotted in figures 3 and 4. The 
remarkable fact is that for the SCG configuration the 
slope is equal to that of the low temperature solutions 
obtained with MFSB, while for the EG there is a small 
renormalization of it. This shows that the essential dif- 
ference between the low and high temperature regime as 
far as the ZS is concerned is the existence of a non-zero 
Y-intercept coordinate at low temperatures as has been 
seen in the experiment. In spite of the similar behavior 



of the slope of the function A(B), the magnetization in 
this high temperature regime is very different. The DMS 
at zero field is much less at high temperature than at low 
temperature. 

Summarizing, we have studied the magnetic response 
of a QD connected with two different geometries SCG 
and EG. To do so we developed a slave-boson formalism 
that takes into account a different boson field for each 
spin. In the Kondo regime there is a critical magnetic 
field B c above which the Kondo resonance can be seen 
splitted. B c depends upon the topology of the system, 
being larger for the EG than for the SCG configuration. 
The Zeeman energy associated to this field is lower that 
the Kondo temperature. For B < B c , A suffers a non- 
linear and rapid growth with the magnetic field that oc- 
curs due to the spin dependent local energy Coulomb 
renormalization and the quenching of the LSF induced 
by the magnetic field. When the spin fluctuations are 
completely quenched, B > B c , the spin splitting is pro- 
portional to the magnetic field with a slope independent 
on temperature and gate potential. The Y-interception 
coordinate is non zero and increases when the gate poten- 
tial is reduced. At a fixed magnetic field, A(V g ) results 
to be a linear function with a B independent slope. The 
value of A is always greater than Ab for T < Tfc due to 
the spins correlations implied in the Kondo effect. Our 
results are in good agreement with recent experimental 
results. 
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